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Abstract 


This  report  documents  a  two-dimensional,  finite  element  study  of  the  evaporation  of  a  liquid  on  a  Hat 
surface  into  a  laminar  air  stream  for  three  geometries:  a  hemispherical  drop  on  a  flat  plate;  a  hemispherical  drop 
between  parallel  plates;  a  smear  on  a  flat  plate.  The  problems  are  made  dimensionless  and  the  flux  from  the 
liquid  is  determined  in  u  parametric  analysis  with  a  range  of  values  of  the  resulting  independent  variables.  A 
comparison  is  made  with  some  experimental  data  and  an  example  of  the  use  of  the  results  is  given. 


Rdsumd 


Cc  rapport  documente  uiiu  dtude  bi-dimensionnelle  utilisant  un  moddle  d’dldmcnt  ddfmi  de  I’dvaporuiion  d’un 
liquidc  sur  une  surface  plane  dans  un  passage  d’air  laminuire  sur  trois  configurations:  Une  guuttu  h(3misphCrk|U0 
sur  une  surface  plane;  Une  gouttc  hdmisphdrique  entre  deux  surfaces  parallilcs;  Une  tache  sur  une  surface  plane. 
Lcs  probldmes  sont  faits  sans  dimensions  et  la  fluctuation  provenant  du  liquide  cst  calculiie  dans  une  unaly.se 
param6trique  ayant  une  dchelle  de  valeurs  pour  les  variables  inddpendantes  rdsultantes.  Une  compuruison  esl 
faite  avec  certaines  donndes  expdrimentales  et  un  exemple  de  I’usage  des  rdsultats  numdriques  cst  donnde. 
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Extoitivc  Summary 


Evaporation  of  liquids  from  fabrics  is  of  interest  in  a  number  of  areas  of  research  such  as  the  study  of 
thermoregulation  by  evaporation  of  sweat  from  the  body,  the  study  of  heat  and  moisture  transport  from  clothing 
and  the  study  of  the  performs  tee  of  protective  clothing  when  challenged  by  drops  of  chemical  agent. 

In  order  to  adequately  predict  the  evaporative  heat  flux  or  the  performance  of  protective  clothing  for 
u  variety  of  conditions,  a  knowledge  of  the  evaporation  rate  of  liquids  from  the  clothing  should  be  known.  This 
information  can  be  obtained  from  a  number  of  sources,  most  of  which  are  complimentary  rather  than  competing 
techniques.  Reid  tests  are  useful  since  they  closely  approidmate  actual  applications,  however,  the  large  number 
of  variables  involved  may  complicate  the  study  to  the  point  where  the  results  are  too  specific  to  extrapolate  to 
other  cases  or  the  effects  cannot  be  attributed  to  individual  causes.  Laboratory  tests  can  be  more  uuslly 
controlled  than  field  tests  and  parameters  may  bo  varied  to  examine  the  effects  of  different  variables  but  there 
remains  the  problem  of  extending  the  laboratory  data  to  field  conditions,  Numerical  analysis  can  be  a  source 
of  detailed  data  which  can  be  substantiated  by  both  laboratory  and  field  data.  Numerical  analyses  cun  provide 
detailed  information  on  individual  aspects  of  a  problem  which  may  not  be  measurable  in  physical  experiments. 
It  can  also  be  used  to  more  confldently  extend  laboratory  data  to  field  conditions.  Thus,  numerical  modelling 
provides  the  investigator  with  a  complimentary  tool  for  understanding  the  physics  and  physical  relationships 
involved  in  a  problem. 

This  paper  reports  a  parametric  study  of  the  transport  characteristics  from  surfaces  into  a  laminar  Utiid 
flow  by  means  of  a  two-dimensional,  flnite  element  analysis  for  three  simple  geometric  configurations.  The 
configurations  which  were  investigated  are  evaporation  of  a  drop  between  two  parallel  plates,  the  evaporation 
of  a  drop  from  a  flat  plate  and  evaporation  from  a  smear  on  a  flat  plate.  Although  the  analysis  was  kept  as 
general  as  possible,  these  geometries  were  selected  as  they  are  encountered  in  the  study  of  protection  from 
hazardous  chemicals.  The  case  of  a  drop  between  two  parallel  plates  is  similar  to  an  experimental  apparatus 
employed  to  measure  the  vapour  penetration  of  a  chemical  agent  drop  on  a  repellent,  chemically-protectivc 
fabric.  The  flat  plate  geometry  Ls  similar  to  the  parallel  plate  geometry,  however,  it  is  more  representative  ol 
a  drop  on  a  repellent,  chemically-protectlve  garment.  The  smear  on  a  flat  plate  geometry  corresponds  to  the 
situation  of  a  drop  on  a  chemlcally-protective  fabric  that  permits  spreading  of  the  liquid  into  the  fabric.  The 
equations  governing  these  problems  are  made  dimensionless,  that  is  without  units,  to  make  the  results  more 
genorally  applicable  to  other,  similar  problems. 

A  parametric  study  using  dimensionless  variables  is  shown  to  provide  a  concise  presentation  of  the 
results  which  can  then  be  used  to  supply  information  on  particular  problems  with  little  additional  work.  Using 
dimensionleu  variables,  both  in  physical  experiments  and  numerical  analyses,  allows  the  investigator  additional 
flexibility  since  the  equations  of  heat  transfer  and  mass  transfer  have  similar  form.  This  technique  also  minimizes 
the  number  of  independent  variables  required  to  characterize  the  problem. 

The  results  presented  in  this  report  can  be  used  to  predict  evaporation  rates  from  drops  (or  heat  fluxes 
from  cylinders)  for  three  simitar  geometries.  The  results  for  the  specific  example  of  a  drop  of  mustard  between 
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parallel  plates  compare  very  favourably  with  measured  values  from  drops  in  experimental  test  cells,  although  thu 
comparison  is  limited  to  one  set  of  experimental  results.  Although  limited,  the  good  agreement  lends 'confidence 
to  the  validity  of  the  numerical  results  and  to  the  use  of  the  results  to  make  predictions  for  cases  difficult  to 
obtain  in  the  laboratory.  Extension  of  the  comparison  to  include  field  data  under  similar  conditions  remains  to 
be  done  and  so  no  conclusion  can  be  made  in  this  area  based  on  the  results  of  this  study. 

An  analysis  of  this  type  would  also  be  of  use  to  designers  of  a  test  apparatus  for  measuring  agent 
penetration  into  protective  fabrics.  By  selecting  the  desired  air  speed  and  drop  size,  calculations  can  be  made 
of  the  appropriate  dlmeniionleu  variables.  The  theoretical  results  can  then  be  used  os  guidelines  for  establishing 
parameters  of  the  apparatus  such  as  channel  height  and  drop  position  so  that  the  apparatus  itself  does  not 
influence  the  outcome  of  the  experiment. 

It  should  bo  stressed  that  the  two-dimensional  study  presented  hero  is  only  an  approximation  to  the 
actual,  three  dimensional  physical  problem  and  therefore  the  results  must  bo  used  with  some  caution,  preferably 
with  experimental  corroboration.  The  comparison  of  a  flow  over  a  cylinder  to  that  over  a  sphere,  which  can  be 
found  in  most  elementary  fluid  mechanics  and  heat  transfer  texts,  further  illustrates  the  caveats  which  must  lie 
observed  when  approximating  a  three-dimensional  problem  with  a  two-dimensional  model.  More  powerful 
computers  would  allow  the  analysis  to  be  extended  into  three  dimensions  which  should  make  the  results  more 
representative  of  the  physical  problem  and  which  could  be  used  to  explore  more  interesting  geometries. 
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1.0  Introduction 


Evaporation  of  liquids  from  fabrics  is  of  interest  in  a  number  of  areas  of  research  such  as  the  study  of 
thermoregulation  by  evaporation  of  sweat  from  the  body,  the  study  of  heat  and  moisture  transport  from  clothing 
und  the  study  of  the  performance  of  protective  clothing  when  challenged  by  drops  of  chemical  agent. 

In  order  to  adequately  predict  the  evnpornlivc  heat  flux  or  the  performance  of  protective  clothing  for 
u  variety  of  condition.^,  a  knowledge  of  the  evaporation  rate  of  liquids  from  the  clothing  .should  be  known.  This 
information  can  be  obtained  from  u  number  of  sources,  most  of  which  are  complimentary  rather  than  competing 
techniques.  Field  tests  arc  useful  since  they  closely  approximate  actual  applications,  however,  the  large  number 
of  variables  involved  may  complicate  the  study  to  the  point  where  the  results  are  too  specific  to  extrapolate  to 
other  coses  or  the  effects  cannot  be  attributed  to  individual  causes.  Laboratory  tests  cun  be  more  ea.sily 
controlled  than  field  tests  and  parameters  may  be  varied  to  examine  the  effects  of  different  variables  but  there 
remains  the  problem  of  extending  the  laboratory  data  to  field  conditions.  Numerical  analysis  can  be  a  source 
of  detailed  data  which  can  be  substantiated  by  both  laboratory  and  Held  data.  Numerical  analyses  cun  provide 
detailed  information  on  individual  aspects  of  a  problem  which  may  not  be  measurable  in  physical  experimunis. 
it  cun  also  be  used  to  more  confidently  extend  laboratory  data  to  field  conditions.  Thus,  numerical  modelling 
provide.s  the  investigator  with  u  complimentary  tool  for  understanding  the  physics  and  phy.sical  relutiomships 
involved  in  a  problem. 

This  paper  reports  a  parametric  study  of  the  transport  characteristics  from  surfaces  into  a  laminar  lluid 
i1ow  by  means  of  a  two'dimensional,  finite  element  analysis  for  three  simple  geometric  configurations,  The 
configurations  which  were  investigated  are  evaporation  of  a  drop  between  two  parallel  plates,  the  evaporation 
of  a  drop  from  a  flat  plate  and  evaporation  from  a  smear  on  a  flat  plate.  Although  the  analysis  was  kept  ns 
general  as  possible,  these  geometries  were  selected  as  they  are  encountered  in  the  study  of  protection  from 
hazardous  chemicals.  The  case  of  a  drop  between  two  parallel  plates  is  similar  to  an  experimental  apparatus 
employed  to  measure  the  vapour  penetration  of  a  chemical  agent  drop  on  a  repellent,  chemically-proteetivc 
fabric.  The  flat  plate  geometry  is  similar  to  the  parallel  plate  geometry,  however,  it  is  more  rcprc.suntutlvc  of 
a  drop  on  a  repellent,  chemically*protectivc  garment.  The  smear  on  a  flat  plate  geometry  corresponds  to  the 
situation  of  a  drop  on  a  chemically-protectivo  fabric  that  permits  spreading  of  the  liquid  into  the  fabric.  I'he 
equations  governing  these  problems  are  made  dimensionless,  that  is  without  units,  to  make  the  results  more 
generally  applicable  to  other,  similar  problems. 
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2.0  Theory 


2.1  Background 

The  equations  governing  the  conservation  of  mas.s,  momentum,  energy  and  chemical  species  ibr  an 
incompressible,  laminar  flow  can  be  written  as  [Currie  1974,  Rohsenuw  1961]; 


Vu“0 

(1) 

(2) 

(3) 

u-Vc^D'^c. 

d  d 

(4) 

Buoyancy  effects  have  been  neglected  and  the  mass  concentration  of  the  chemical  species  is  assumed  to  be 
sufficiently  small  that  it  docs  not  significantly  affect  the  fluid  properties.  Definition  of  the  terms  in  the  equations 
cun  be  found  in  the  glossary.  These  equations  can  be  non-'dimensionalized  through  the  introduction  of 
appropriate  geometric  and  fluid  property  constants  relevant  to  the  problem.  These  constants  typically  inuliKlu: 
a  length  scale  such  as  a  boundary  length  or  channel  width;  a  velocity  scale  such  as  the  mean  inlet  velocity;  a 
temperature  scale  such  as  the  difference  between  the  inlet  fluid  temperature  and  a  hot  wall  temperature;  and 
a  chemical  species  mass  concentration  scale  such  as  the  difference  between  the  chemical  species  concentration 
of  the  inlet  fluid  and  a  source  on  another  boundary.  These  constants  in  addition  to  physical  properties  of  the 
lluid  such  as  vi.scosity,  thermal  conductivity  and  mu.ss  diffusivity  can  be  used  to  group  a  large  number  of  variables 
into  a  much  smaller  number  of  non-dimensional  variables  which  characterize  the  flow,  This  greatly  simplillcs 
a  parametric  study  of  a  given  problem.  Proper  non-dimensionalization  of  the  governing  equations  also  gives  the 
investigator  an  indication  of  which  terms  are  significant  in  the  analysis  once  numerical  values  are  assigned  to  the 
various  constants  and  fluid  properties. 

In  dimensionless  form,  the  equations  may  be  written  as: 

V‘l/-0  (5) 


(6) 

Re 


i/.vr«-^v^  (7) 

RePr 
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I/.VC=--i-X^C’ 

ReSc 


(8) 


Note  the  similarity  between  the  energy  transport  equation  and  the  mass  transport  equation.  By  simply  changing 
the  non'dimensional  quantity  from  RePr  to  ReSc,  results  from  a  heat  transfer  experiment  (either  numerical  or 
physical)  can  be  used  to  describe  mass  transfer  in  the  same  geometry.  This  can  be  useful  as  heat  transfer 
experiments  arc  frequently  easier  to  perform  than  ma.ss  transport  experiments  and  yet  will  provide  essentially 
equivalent  information. 


The  non-dimensional  heat  or  mass  transfer,  Q,  from  a  surface  can  be  converted  back  into  a  dimensional 
quantity,  q,  through  the  following  equations; 


9  = 


or 


(10) 


respectively, 

The  boundary  conditions  required  to  complete  the  mathematical  statement  of  a  problem  are  the 
specificution  of  velocity,  temperature  and  chemical  species  along  all  of  the  domain  boundaries  or  u  coinbinution 
of  these  primitive  quantities  along  some  of  the  domain  boundaries  and  specification  of  fluxes  along  the  remaining 
boundaries.  In  a  non-dimensional  analysis,  appropriate  dimensionless  boundary  conditions  must  be  supplied. 

This  study  will  only  consider  a  laminar,  isothermal,  mass  transport  problems  using  equations  5,  6,  S  and 
10  with  terminology  consistent  with  mass  transport  problems.  The  results  are,  as  noted  before,  equally  applicable 
to  heat  transfer  problems  of  similar  geometries  without  a  chemical  contaminant  present.  It  should  be  noted  that, 
when  considering  a  heat  transfer  problem,  radiant  heat  transfer  must  also  be  considered  as  a  boundary  condition 
if  the  flowing  fluid  is  a  gas. 


2JI  Problem  Description 

The  three  geometries  studied  are  shown  in  Figure  2.2.1  with  an  accompanying  description  of  the 
boundary  conditions  and  the  associated  non-dimensionalization  parameters  used  for  each  analysis. 

Figure  2.2.1a  shows  a  drop  at  some  position  downstream  of  the  leading  edge  of  a  flat  plate.  Variable.s 
of  importance  in  this  problem  are  the  drop  radius,  r,  the  boundary  layer  thickness  at  the  drop,  S,  the  mean  How 
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X - - 

— 1 


_ (C) _ I 

Figure  2.2.1  Problem  geometries  for  the  study  of  mass  transfer  in  laminar  air  flows;  (a)  Flow  over  a  drop 
on  a  flat  plate;  (b)  Flow  over  a  drop  between  parallel  plates;  (c)  Flow  over  a  smear  on  a  flat 
plate. 


velocity,  Uq,  the  diffusivity  of  the  drop  species  in  the  fluid,  D,  the  fluid  viscosity,  p,  and  the  fluid  density,  p.  Thu 
boundary  layer  thickness  is  a  function  of  the  distance  along  the  plate,  x,  and  so  the  boundary  layer  thickness  ut 
the  drop  may  be  replaced  by  the  distance  from  the  leading  edge  of  the  plate  to  the  drop.  All  of  thesu  , 
parameters  can  be  combined  into  three  dimensionless  groups  which  can  be  used  to  characterize  the  flow; 


X  -  x/r 

(H) 

Re  - 

(12) 

Sc  ■  m/pD 

(13) 

the  dimensionless  drop  position, 
the  flow  Reynolds  Number, 
the  flow  Schmidt  Number, 
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Figure  2.2.1b  shows  u  similar  problem  where  the  drop  is  now  between  two  paraliel  plates.  Thi.s 
introduces  an  additional  parameter  to  the  previous  list,  the  distance  between  the  plates,  h,  which  can  be  non- 
dimensionalized  by  the  drop  radius: 

the  dimensionless  plate  separation,  H  =  h/r  (14) 

Also,  the  boundary  layer  growth  is  now  limited  to  one-half  the  distance  between  the  two  plates  except  in  the 
vicinity  of  the  drop  which  disturbs  the  How.  If  the  drop  is  within  the  fully  developed  portion  of  the  flow  (where 
the  boundary  layer  extends  to  the  middle  of  the  plates),  the  results  should  not  change  as  the  drop  position  along 
the  plates  is  increased.  Thus,  the  drop  position  along  the  plate  is  only  important  in  the  initial,  developing  region 
of  the  flow. 

Figure  2.2.1c  shows  a  somewhat  simpler  problem,  that  of  a  smear  on  a  flat  plate.  Now,  the  flow  is  not 
dLsturbed  by  the  size  of  the  "drop"  although  the  evaporation  rate  (per  unit  area)  will  depend  somewhat  on  the 
smear  radius.  In  this  problem,  the  plate  length  is  the  more  relevant  length  scale  and  the  problem  parameters 


the  dimensionless  smear  position. 

X  =  x/t 

(15) 

the  dimensionless  smear  size. 

R  -  r/t 

(16) 

the  flow  Reynolds  Number, 

Re  • 

(17) 

the  flow  Schmidt  Number, 

Sc  =■  li/pD 

(1«) 

The  Schmidt  Number  Is  an  indicator  of  the  relative  importance  of  the  diffusivity  of  momentum  to  the 
chemical  diffusivity  in  the  flow,  similar  to  the  Prandtl  number  In  heat  transfer  problems.  In  this  study,  three 
Schmidt  numbers  were  used:  0.6,  0.7  and  2.5.  The  first  two  Schmidt  Numbers  produce  virtually  the  sumo  results 
and  cun  be  considered  equivalent  for  most  purposes.  These  Schmidt  Numbers  arc  representative  of  the  diffusion 
of  water  vapour  in  air  (0.6)  and  the  diffusion  of  mustard  vapour  in  air  (2.5).  A  Schmidt  Number  of  0.7  produces 
results  equivalent  to  those  obtained  from  a  heat  transfer  analysis  into  the  air  stream  with  the  same  geometry 
(that  is  with  a  Prandtl  Number  of  0.7). 

Although  the  emphasis  of  this  study  is  on  the  evaporation  of  a  liquid  into  the  adjacent  laminar  air- 
stream,  the  diffusion  of  the  liquid  contaminant  from  the  drop  through  the  fabric  upon  which  it  rests  is  also  of 
interest.  A  limited  investigation  of  this  is  presented  here  using  a  constant  value  of  the  mass  diffusivity  through 
the  fabric.  In  this  case,  the  fabric  thickness,  f,  is  the  principle  variable. 

With  the  above  transformation,  the  problem  is  to  quantify  the  non-dimensional  ma.ss  flux,  Q,  in  terms 
of  the  relevant  non-dimensional  quantities: 

Q  -  f(X,Rc,Sc), 

Q  -  f  (X,  H,  Re,  Sc), 

Q  -  (^(X,  R,  Re,  Sc), 

Q  -  f"(F), 
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for  a  drop  on  a  plate; 

(ly) 

for  a  drop  between  parallel  plates; 

(20) 

for  a  smear  on  a  plate; 

(21) 

for  diffusion  through  a  fabric. 

(22) 

I 


This  is  accomplished  by  holding  all  but  one  independent  variable  in  the  above  relationships  constant  and 
examining  the  elTcct  of  changing  the  remaining  variable  on  the  flux.  Each  independent  variable  i.s  examined  over 
a  range  which  includes  practical  limits  so  that  the  results  can  be  interpolated  to  yield  information  on  specii'ic 
configurations. 


2.3  Solution  Technique 

Unfortunately,  no  general  solution  exists  for  equations  5  through  8  so  that  it  is  necessary  to  make  use 
of  an  approximation  technique  to  obtain  an  estimate  of  the  solution.  The  technique  used  in  this  study  was  that 
of  finite  elements,  implemented  through  commercial  software  called  Pldap  [Engleman  1991], 

In  the  finite  clement  technique,  the  problem  domain  is  divided  into  many  smaller  domains  or  finite 
elements.  The  governing  equations  are  then  approximated  over  each  of  these  finite  elements.  The  solution  at 
the  various  nodes  within  these  elements  is  obtained  by  specifying  the  form  of  the  di.strlbution  of  the  independent 
variables  over  the  elements,  re-assembling  the  elements  into  a  global  matrix  and  then  solving  the  resulting 
simultaneous  nonlinear,  algebraic  equations.  In  this  study,  nine-node  quadrilateral  elements  were  used. 

Two  different  geometries  were  considered  for  the  solid  fabric  surface  below  the  liquid  contaminant.  In 
the  two  problems  involving  a  drop  (Figure  2.3.1a),  it  was  assumed  that  the  drop  rested  on  a  liquid  repellent  but 
vapour  permeable  surface.  Thus,  diffusion  of  the  vapour  through  the  lower  boundary  was  included  in  the 
analysis.  In  the  remaining  problem,  it  was  assumed  that  the  boundary  was  not  repellent  but  instead  promuied 
spreading  of  the  liquid  resulting  in  a  smear.  In  this  problem,  diffusion  within  the  porous  surface  wius  not 
modeled  for  simplicity  (Figure  3,2.1b).  It  should  be  noted  that  lateral  diffusion  in  the  fabric  would  probably  be 
present  in  an  actual  application. 

Because  of  the  non-linear  nature  of  the  governing  equations,  the  solution  is  found  using  an  iterative 
procedure.  Typically,  the  iterative  procedure  involved  between  two  to  six  steps  of  successive  substitution  followed 
by  as  many  as  forty  steps  of  Quasi-Newton  iterations.  Acceleration  factors  between  0.5  and  0.9  were  used  fur 
problems  in  which  convergence  was  difficulc.  Streamwise  Upwinding,  a  technique  frequently  used  to  promote 
stability  of  the  numerical  analysis  was  not  used  in  the  solution.  Convergence  was  assumed  when  each  of  the 
following  criteria  were  met:  the  norm  of  the  relative  change  of  the  solution  vector,  |U|  •  Uj.il/|U||,  was  less 
than  0.001;  the  norm  of  the  relative  change  of  the  residual  vector,  IRjI/IRqI  was  less  than  O.iXll;  the  maximum 
norm  of  the  residue,  |U{-Uol,  was  less  than  approximately  0.05. 


3.0  Results  and  Discussion 
3.1  Drop  on  a  Plate 

Figure  3.1.1  shows  u  typical  rc.sult  for  the  flow  field  streamlines  and  mass  concentration  contours  i'or 
transport  from  a  drop  on  a  fiat  plate  in  a  laminar  flow  field.  Most  of  the  evaporated  mass  remains  close  to  the 
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Figure  2.3.1  Sections  of  the  finite  element  meshes  for  flow  analyses:  (a)  flow  over  a  drop  on  a  flat  plate 


and  a  drop  between  parallel  plates;  (b)  flow  over  a  smear  on  a  fiat  plate. 

plate;  little  is  transported  far  into  the  flow  which  is  typical  of  laminar  transport  flows.  The  vapour  from  the  drop 
can  also  be  seen  to  spread  laterally  through  the  fabric  for  approximately  four  drop  radii  downstream  and  three 
drop  radii  upstream  of  the  drop  edges. 

Figure  3.1.2  shows  the  resulting  mass  flux  into  the  air-stream  as  a  function  of  position  along  the  plate 
for  two  Schmidt  Numbers  at  a  selected  drop  Reynolds  Number  of  30.  The  flux  decreases  with  increasing 
distance  from  the  leading  edge  of  the  plate,  or  as  the  boundary  layer  thickness  increases.  The  shape  of  the 
curve  suggests  an  asymptotic  decrease  with  increasing  distance.  For  a  drop  Reynolds  Number  of  30,  little  change 
in  the  flux  occurs  beyond  L=50.  The  position  along  the  plate  after  which  the  flux  remains  effectively  constant 
will  probably  increase  somewhat  with  an  increasing  drop  Reynolds  Number,  although  this  is  speculation  as  this 
study  did  not  examine  this  point. 

Figure  3.1.3  shows  the  mass  flux  from  the  drop  into  the  air-stream  as  a  function  of  the  drop  Reynolds 
Number  at  two  positions  along  the  plate  and  for  two  Schmidt  Numbers.  As  the  drop  Reynolds  Number 
increases,  the  mass  flux  also  inaeases  although  the  changes  are  diminishing  as  the  drop  Reynoids  Number 
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Figure  3.1.1 


Typical  slrcumlincs  (.solid  lines)  and  concentration  contours  (dashed  lines)  for  How  over 
a  drop  (solid  hemisphere)  on  a  flat  plate. 
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Figure  3.1.2  Flux  from  the  upper  surface  of  a  drop  into  the  airstream  as  a  function  of  position  along  the 
plate  for  two  different  Schmidt  Numbers  and  a  drop  Reynolds  Number  of  30. 


increases.  The  flux  will  likely  increase  abruptly  and  substantially,  however,  when  the  flow  becomes  turbulent  as 
the  rate  of  convective  mass  transport  increases  significantly  with  the  onset  of  turbulence. 

It  should  be  noted  that  an  increasing  drop  Reynolds  Number  can  result  either  from  an  increasing  drop 
size  or  an  increasing  flow  velocity  and  that  the  dimensional  mass  flux  is  inversely  proportional  to  the  drop  size 
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Figure  3.1.3  Flux  from  the  upper  .surface  of  a  drop  on  a  flat  plate  as  u  function  of  drop  Roynold.s 
Number  for  two  positions  along  the  plate  and  two  .Schmidt  Numbers. 


for  thi.s  probl  -mi 

(1  -  Q(DAc„/r)  (2.3) 

Thu.s,  even  though  the  non-dimonsionul  flux  may  increase  with  an  increasing  drop  Reynolds  Number,  the 
dimensional  mass  flux  may  be  quite  different  depending  on  whether  the  larger  Reynolds  Number  change  is  due 
to  increasing  drop  si/«  or  increasing  flow  velocity. 

The  cases  with  ihc  large  Schmidt  Number  have  the  greater  non-dimensional  mass  flux  and  again  the 
mass  flux  is  greater  closer  to  the  leading  edge  of  the  plate  where  the  boundary  layer  is  thinner.  At  Hrst  glance 
this  is  somewhat  surprising  as  the  Schmidt  Number  is  inversely  proportional  to  the  chemical  diffusiviiy,  however, 
on  inspection,  tiie  dependence  of  the  non-dimensional  flux  is  not  linearly  proportional  to  the  Schmidt  Number 
and  when  actual  mass  flaxes  are  calculated,  using  equation  23,  the  chemical  species  with  the  greater  diffusivlty 
does  indeed  have  the  greater  flux. 

As  expected,  the  mass  flux  downwards  (through  the  fabric)  showed  no  dependence  on  the  position  along 
the  plate.  As  only  one  value  of  the  non-dimensional  plate  thickness  was  u.sed  in  this  portion  of  the  investigation, 
no  general  conclusions  can  be  drawn.  It  is  expected  that  the  downwards  flux  for  this  geometry  will  be  virtually 
identical  to  that  for  a  drop  in  a  channel  which  is  discussed  in  the  next  section, 
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3.2  Drop  In  A  Channel 


Figure  3.2.1  .shows  lypical  streamlines  and  concentration  contours  for  flow  over  a  drop  between  paruilel 
plute.s.  Notice  that  there  is  some  horizontal  diffusion  from  the  bottom  of  the  drop  within  the  lower  plate  which 
then  diffuses  upwards  into  the  uir.strcum  both  in  front  of  and  behind  the  drop. 


Figure  3.2.1  Typical  .struumllnes  (solid  lines)  and  concentration  contours  (dashed  lines)  for  tlow  over  a 
drop  (solid  hemisphere)  in  a  channel. 


The  vertical  flux  showed  only  a  slight  dependence  on  the  position  along  the  plates  (Figure  3.2.2)  with 
the  flux  being  marginally  greater  near  the  leading  edge  of  the  channel.  The  results  indicate  that  if  the  non- 
dimensional  distance,  L/r,  from  the  leading  edge  of  the  channel  is  greater  than  approximately  10,  no  signiricanl 
dependence  on  position  will  be  observed  for  drop  Reynolds  Numbers  less  than  30. 

Only  a  very  slight  dependence  of  the  vertical  flux  on  the  channel  height  was  found  (Figure  3.2.4).  As 
the  drop  height  approaches  that  of  the  channel  height,  the  flow  becomes  disturbed  by  the  drop,  especially  in  a 
two-dimensional  analysis,  and  the  flux  increases  accordingly.  The  results  indicate  that  if  the  height  of  the  channel 
is  greater  than  10  to  15  times  that  of  the  drop  size,  then  the  channel  height  will  have  a  negligible  influence  on 
the  vertical  flux  from  the  drop. 

The  vertical  flux  increases  with  an  increasing  drop  Reynolds  Number.  Figure  3.2.3  show^  this 
dependence  at  two  positions  along  the  plate,  5  and  50. 

Figure  3.2.5  shows  the  downwards  mass  flux  due  to  diffusion  through  the  porous  (fabric)  wall  as  u 
function  of  the  wall  thickness.  As  one  would  expect,  all  flow  parameters  (Re,  H/r,  L/r)  other  than  the  wall 
thickness  have  no  effect  on  the  downward  flux  which  b  shown  in  Figure  3.2.5  where  all  data  fall  on  a  single 
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Figure  3.2.2  Flux  from  the  upper  surface  of  u  drop  as  a  function  of  position  along  the  channel  for  two 
Schmidt  Numbers.  The  Drop  Reynolds  Number  was  30  and  the  channel  height  was  .V 


Plate  Reynolds  Number,  Re 


Figure  3.2.3  Rux  from  the  upper  surface  of  a  drop  in  a  channel  as  a  function  of  drop  Reynolds 
Number  for  two  Schmidt  Numbers  and  two  positions  along  the  channel.  The  channel 
height  was  5. 


curve.  This  also  implies  that  the  downwards  flux  from  a  drop  on  a  flat  plate,  described  in  the  previous  section, 
should  be  identical  to  the  downwards  flux  in  a  channel  for  the  same  wail  thickness  and  concentration  difference. 
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Figure  3.2.4  Flux  from  ihu  upper  surface  of  a  drop  In  a  chunnol  as  a  function  of  the  channel  height. 

Two  Schmidt  numbers  and  two  positions  are  shown.  The  drop  Reynolds  Number  wus  30. 


Figure  3.2.5  The  flux  from  the  lower  surface  of  a  drop  in  a  channel  as  a  function  of  the  wall  thickness. 


The  Insert  gives  the  results  In  greater  detail  for  thin  walls. 

As  the  wall  thickness  increases,  the  flux  decreases  due  to  the  decreasing  concentration  gradient  across 
the  wall.  The  results  suggest  that  as  the  wall  thickness  decreases  to  zero,  the  mass  flux  becomes  very  large  If 
not  Infinite  and  as  the  wall  thickness  increases,  the  flux  tends  to  zero.  The  greatest  change  in  the  downwards 
mass  flux  occurs  in  the  wall  thickness  range  of  zero  to  fifteen.  For  wall  thicknesses  greater  than  fifteen,  a 
substantial  increase  in  the  wall  thickness  is  required  in  order  to  noticeably  reduce  the  mass  flux. 
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33  Smeur  On  A  Flat  Piute 


Figure  3.3.1  shows  a  typical  set  of  streamlines  and  concentration  contours  for  the  How  over  a  smeur  on 
u  ilut  plate.  As  diffusion  within  the  plate  was  nut  modeled  in  this  problem,  there  is  no  mass  flux  into  the  Huid 
iK'fore  the  start  of  the  drop.  Bu.sed  on  the  results  of  the  previous  section,  one  would  expect  .some  dilTu.sion  and 
muss  flux  from  the  wall  into  the  Huid  to  occur  both  before  and  after  the  contaminulcd  region, 


Figure  3.3.1  Typical  streamlines  (solid  lines)  and  concentration  contours  (du.shcd  lines)  for  lluw  over 
u  .smeur  (thick  solid  line)  on  u  flat  plate. 


The  vertical  mass  flux  from  the  smear  was  found  to  vary  with  position  along  the  plate  (Figures  3.3.2  and 
3.3.3),  although,  most  of  the  differences  occur  within  the  first  quarter  of  the  plate  length.  The  variation  is  most 
pronounced  at  the  lower  Reynolds  Numbers  as  the  boundary  layer  increases  in  thickness  more  rapidly  in  these 
cases  and  the  convective  mass  transport  is  much  less  intense. 

The  vertical  moss  flux  varies  considerably  with  an  increasing  drop  Reynolds  Number,  especially  beyond 
the  first  quarter  of  (he.  plate  length  where  the  boundary  layer  thickness  is  changing  as  the  inverse  .squure-rool 
of  the  plate  Reynolds  Number.  These  variations  arc  primorily  a  result  of  changes  in  the  fluid  velocity. 


3.4  Use  of  Results 

In  order  to  demonstrate  how  the  results  of  dimensionless  studies  may  be  used,  an  example  is  presented 
for  each  of  the  geometries  analyzed  above.  It  is  as.sumcd  that  a  1  microlitre  drop  of  liquid  mustard  (HD)  having 
u  diffusivity,  D,  of  6x10*^  m^/s  is  placed  at  a  position,  (,  2.5  cm  from  the  leading  edge  of  a  5  cm  long  plate.  The 
drop  is  exposed  to  a  mean  air  speed,  u,  of  0.5  m/s  and  the  air  temperature  is  30*C.  For  the  flow  between 
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Reynolds  Numbers.  The  Schmidt  Number  was  2.5, 

parallel  plates  geometry,  the  plate  separation,  h,  is  0.5  cm.  For  the  case  of  a  smeur  on  a  flat  plate,  the  fabric 
is  assumed  to  be  0,5  mm  tliick  and  the  volume  of  fabric  through  which  the  drop  spreads  is  half-filled  with  fabric 
threads.  The  volatility  of  the  liquid  is  assumed  to  be  approximately  1.2  g/m^  and  the  inlet  air  is  assumed  to  have 
a  2cro  concentration  of  the  mustaid  vapour. 
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Figure  3.3.4  Typical  (rend.i  of  the  flux  from  a  smear  on  a  flat  plate  us  a  function  of  pu.sitiun  at  two 
Schmidt  Numbers.  The  Drop  Reynolds  Number  was  2866. 


Noti-dimensiunulizing  the  physical  parameters  for  a  drop  on  a  flat  plate  yields  the  following; 

2n 


Sr-  ,73^75 

(1.3itg/m’)(fixiO-^wVi) 

^  _  (1.3i^/m^K0.5i>t/5)(7.8x70-*w) 
i.SxlO~*kglm»s 


(25) 


(26) 


£^:.2;5£i0jm_j2  (27) 

7.9xl0'*m 

From  figure  3.1.4,  the  dimensionless  flux  into  the  airstream  from  the  upper  surface  of  the  drop  is  found  to  be 
approximately  8.5.  This  can  be  converted  back  mto  a  dimensional  flux: 
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.omglmh 

l.SxW*m 


(28) 


which,  when  multiplied  by  the  upper  surface  area  uf  the  drop  ^ves  un  absolute  flux  of  3x10*^  g/s 


For  the  same  drop  between  parallel  plates,  most  of  the  above  dimensionless  variables  can  be  used  along 
with  the  dimensionless  plate  separation; 


„  O.SxIO'^m  ^ 

fjm - - —  -5 

lAxlO'*m 


(29) 


From  Figure  3.2.3,  the  dlmen.sionle.ss  flux  is  found  to  be  approximately  9.  This  is  equivalent  to  u  dimensionul 
llux  of  0.083  g/m^s  or  un  absolute  flux  of  3, 2x10'’  g/s.  This  Is  similar  to  that  found  for  the  llux  from  a  drop  on 
u  Hut  plate  which  is  to  be  expected  us  H,  the  ratio  of  the  channel  height  to  the  drop  height,  is  large. 

If  the  liquid  soaks  into  the  fabric,  the  resulting  smear  radius  is: 

r « 1.  UiO -»m 
(0.5jei0’^m)ii 

The  plate  Reynolds  Number  is; 

(1.3<;g/m^)(0.5m/g)(0.05m) 
l.SxiO  ■*%/»!•« 

The  dimensionless  position  is  H-0,S.  From  Figure  3.3.3,  the  dimensionless  flux  is  found  to  be  approximately 
2  which  can  be  converted  into  dimensionai  form  as: 

q  .2  ~0.OQO3glm  (32) 

0.05m 

The  spread  uf  the  drop  on  and  through  the  fabric  depends  a  great  deal  on  the  fabric  inatcriul  and  on  the 
construction  both  of  the  fabric  itself  (its  weave)  and  of  the  individual  yarns  of  the  fabric  [Crow  1989].  If 
extrapolation  from  experimental  data  [Crow  1989]  to  smaller  drop  volumes  is  valid,  a  1  pi  drop  should  spread 
between  1x10*^  suid  4x10'^  which  yields  a  very  small  absolute  flux  of  10'^  g/s.  If  the  smear  spreads  along  the 
yarns  of  the  fabric  inster  i  of  just  fllling  the  voids  between  the  yarns,  the  area  of  the  drop  could  be  a  good  deal 
larger  and  the  resulting  flux  more  substantial,  however,  the  smear  would  have  to  increase  its  area  two  orders  of 
magnitude  to  be  comparable  to  the  evaporation  from  the  dropa. 


(30) 


(31) 
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3.5  Comparison  of  Numerical  and  Experimental  Results 


Measurements  of  the  evaporation  of  mu.<(tard  (HD)  drops  have  been  made  in  u  standard  Penetration 
Test  Cell  [Fagotto,  1901].  Five,  onc-microlltrc  drops  were  distributed  over  a  treated  fabric  and  exposed  to  un 
air  flow  velocity  of  approximateiy  0.44  m/s  ut  30°C.  Due  to  the  repeilent  treatment  on  the  fabric,  the  drop.s 
remained  approximately  hemispherical.  Thc.se  conditions  are  similar  to  the  example  conditions  for  flow  over  u 
drop  between  parallel  plates  discussed  in  the  previous  section. 

Air  contaminated  with  vapour  leaving  the  cell  was  sampled  by  a  gas  chromatograph  and  then  run 
through  a  bubble  extractor,  The  instantaneous  concentration  of  vapour  in  the  air  was  measured  once  every 
minute.  From  this,  an  average,  instantaneous  evaporation  rate  was  calculated.  As  a  check  on  the  accuracy,  these 
instantaneous  evoporation  rates  were  time  integrated  to  determine  the  total  mass  evaporated  from  the  drops. 
This  result  was  then  compared  to  the  total  measured  by  the  bubble  extractor  and  the  results  where  found  to  be 
in  agreement  to  within  5%.  Since  the  bubble  extractor  measures  only  the  cumulative  total  ma!i.s  of  vapour 
leaving  the  cell,  it  does  not  directly  corroborate  the  instantaneous  evaporation  rate,  however,  it  does  lend 
credence  to  the  results. 

After  un  initial  period  during  which  the  apparatus  came  up  to  temperature,  the  average,  Instuntuneous 
cvuparutlon  rate  per  drop  seemed  to  achieve  u  rclnllvcly  constant  2.7x10'^  g/s.  This  rate  continued  until 
approximately  8t)%  of  total  liquid  hud  evaporated  after  which  the  evaporation  rate  dccrea.sed  steadily, 

The  measured  evaporation  rate  (2, 7x10'^  g/s)  compares  qullc  favourably  to  the  predicted  rate  I'lniiul  in 
the  previous  section  (2, 4x10'^  g/s)  which  was  approximately  12%  Ic-ss  than  the  measured  value. 


4.0  Conclusion 

A  parametric  study  using  dimensionless  variables  Is  shown  to  provide  a  concise  presentation  of  the 
results  which  can  then  be  used  to  supply  information  on  particular  problems  with  little  additional  work.  Using 
dimensionless  variables,  both  in  physical  experiments  and  numerical  analyses,  allows  the  investigator  additional 
flexibility  since  the  equations  of  heat  transfer  and  mass  transfer  have  similar  form.  This  technique  also  minimizes 
the  number  of  independent  variables  required  to  characterize  the  problem. 

The  results  presented  in  this  report  can  be  used  to  predict  evaporation  rates  from  drops  (or  heat  fluxes 
from  cylinders)  for  three  similar  geometries,  The  results  for  the  specinc  example  of  u  drop  of  mustard  between 
parallel  plates  compare  very  favourably  with  measured  values  from  drops  in  experimental  test  cells,  although  the 
comparison  is  limited  to  one  sot  of  experimental  results.  Although  limited,  the  good  agreement  lends  confidence 
to  the  validity  of  the  numerical  results  and  (o  the  use  of  the  results  to  make  predictions  for  cases  difficult  to 
obtain  In  the  laboratory.  Extension  of  the  comparison  to  include  Held  data  under  similar  conditions  remains  to 
be  done  and  so  no  conclusion  can  bo  made  In  this  area  based  on  the  results  of  this  study. 

An  analysis  of  this  type  would  also  be  of  use  to  designers  of  a  test  apparatus  for  measuring  agent 
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penetration  into  protective  fabrics,  By  selecting  the  desired  air  speed  and  drop  size,  calculations  can  be  made 
of  the  appropriate  dimensionless  variables.  The  theoretical  results  can  then  be  used  as  guidelines  for  establishing 
parameters  of  the  apparatus  such  as  channel  height  and  drop  position  so  that  the  apparatus  itself  does  not 
influence  the  outcome  of  the  experiment. 

It  should  be  stressed  that  the  two-dimensional  study  presented  here  is  only  an  approximation  to  the 
actual,  three  dimensional  physical  problem  and  therefore  the  results  must  be  used  with  some  caution,  preferably 
with  experimental  corroboration.  The  comparison  of  a  flow  over  a  cylinder  to  that  over  a  sphere,  which  cun  be 
found  in  most  elementary  fluid  mechanics  and  heat  transfer  texts,  further  illustrates  the  caveats  which  must  be 
observed  when  approximating  a  three-dimensional  problem  with  a  two-dimensional  model.  More  powerful 
computers  would  allow  the  analysis  to  be  extended  into  three  dimensions  which  should  make  the  results  more 
representative  of  the  physical  problem  and  which  could  be  used  to  explore  more  interesting  geometries. 
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Gloiiaury 


Cj,  chcmicul  cuncuntrution 

C,  dimensionless  chemicui  concentration,  c,/ACj 

Ac,  constant  chemical  concentration  difference 
Cp  heat  capacity 

D  thermal  diffusivity 

f  fabric  thickness 

F  dimensionless  fabric  thickness,  f/L 

h  plate  separation 

H  dimensionless  plate  separation,  h/r 

k  thermal  conductivity 

(  plate  length 

L  generic  length  scale  or  dimensionless  plate  length 

p  fluid  pressure 

P  dimensionless  fluid  pressure,  p/C/jpu^^) 

Pr  Prandtl  Number,  p^p/k 

q  heat  or  mass  flux 

0  dimensionless  heat  (qL/kAt)  or  mass  (qL/DAe,)  llux 
r  drop  radius 

R  dlmonsionlesii  drop  radius,  r/L 

Re  Reynolds  Number,  pu^L/g 

Sc  Schmidt  Number,  nlpD 

t  temperature 

T  dimensionless  temperature,  t/Al 

At  constant  temperature  difference 

u  velocity  vector 

U  dimensionless  velocity  vector,  u/u^ 

Ug  reference  speed 

X  distance  from  the  leading  edge  of  a  plate 
X  dimensionless  distance  from  the  leading  edge  of  a  plate,  X/L 

{  boundary  layer  thickness 

fi  fluid  density 

g  fluid  viscosity 

V  gradient  operator 

Laplacian  operator 
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